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We have considered the effects of non magnetic impurities and interface roughness on the interlayer 
coupling between magnetic layers in metallic multilayers. The two types of defects alter the interlayer 
coupling in quite different ways. Elastic electron scattering by impurities in the non magnetic spacer 
layers between magnetic layers produces an exponential decay of the coupling with a characteristic 
decay length that is considerably longer than the "global" transport mean free path for the spacer 
layer with its surrounding interfaces. Interfacial roughness leads to an attenuation of the coupling 
that is related to the width of the roughness in relation to the Fermi wavelength; roughness does not 
alter the range dependence of the coupling. For certain types of electrical transport, e.g., for current 
perpendicular to the plane of the layers, the scattering from interface roughness and impurities in 
the spacer layers contribute on an equal footing to the exponential decay of the electron propagators, 
i.e., global mean free path. We show that interface roughness and impurities in the spacer layer 
affect the interlayer coupling differently. 
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I. INTRODUCTION 

The range dependence of the indirect coupling between 
magnetic ions mediated by conduction electrons, i.e., the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction, in 
the presence of scattering by non magnetic impurities has 
been debated over the past 35 yearsJjia Recently this 
question has been resolved; it is now clear that scatter- 
ing by non magnetic impurities does not alter the range 
dependence of the couphng between two magnetic ions; 
it only introduces a phase factor that depends on the SB& 
cific distribution of impurities between the two ions.uEl 
As long as one does not average over different impurity 
configurations the coupling between a pair of magnetic 
ions remains undamped. Only when one looks at proper- 
ties that require an average of the coupling between two 
ions for different realizations of impurity distributions, 
and therefore by taking an average over the phase factor 
produced by the impurity scattering, does one find that 
the coupling is damped Jd 

The recent interest in the interlayer exchange coupling 
in magnetic multilayersQ raises the questions whether 
scattering in the non magnetic spacer layers dampens 
the coupling between magnetic layers, and what effect 
the roughness of interfaces will have on the range of the 
coupling. Here we show that for the coupling between 
two sheets of spins one does average over different realiza- 
tions of impurity distributions in the intervening spacer. 



and the coupling is exponentially damped as a function 
of the distance between planes, i.e., magnetic layers. The 
damping is proportional to the strength of the impurity 
scattering in the spacer; it has nothing to do with the 
scattering due to roughness at the interfaces. While this 
characteristic decay distance of the interlayer coupling 
and the transport mean free path in the spacer layerQ 
are both due to impurity scattering, they are in no way 
simply related to one|-another. Interfacial roughness at- 
tenuates the coupling^ in proportion to the size of the 
inter diffused region relative to the Fermi wavelength of 
the conduction electrons providing the coupling. How- 
ever, this decrease does not alter the range dependence 
of the coupling. For this reason it would be completely 
erroneous to combine the effects of the scattering from 
impurities in the spacer and interface roughness into one 
decay coefficient to produce an exponential decay of the 
interlayer coupling. While this procedure is correct for 
certain types of electrical transport in magnetic multi- 
layers, e.g., fot; current perpendicular to the plane of the 
layers (CPP),cl interlayer coupling is not in this class of 
situations. Therefore the range of the interlayer coupling 
is considerably longer (characteristic decay of the cou- 
pling is considerably slower) than what one anticipates 
from the transport properties (resistivity) of a magnetic 
multilayer. 

In the following section we first derive the RKKY cou- 
pling between a pair of magnetic ions in the presence of 
non magnetic impurities. Next we show how averaging 
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over the phase induced by the impurity scattering pro- 
duces a couphng that decays with the distance between 
the ions, and the extent to which this can be described by 
an exponential of the distance. We use the approxima- 
tion of representing the coupling between two magnetic 
layers as that due to the two planes of magnetic, ions 
that interface with the nonmagnetic spacer layer,ll3 and 
we show that when one averages over the magnetic ions 
in these planes that this coupling has a slowly decaying 
exponential component as a function of the thickness of 
the spacer layer. Finally we take into account of effect of 
the roughness of the interfaces and show how this atten- 
uates the coupling. 



II. TWO MAGNETIC IONS 

To calculate the coupling between two magnetic ions 
when the conduction electrons are scattered by impu- 
rities we use the approach adopted by Bulaevskii and 
Panyukov.0 It consists of using a semiclassical form for 
the electron propagator (Green's function) connecting 
the positions of the two ions, and taking account of the 
scattering by impurities through the phase of the prop- 
agator. The coupling between two magnetic impurities 
located, respectively, at r and r' is given by 



J(r,r') 



dzTT[G{r,r';z)Gir',r;z)]. (1) 



Here, G is the Green's function corresponding to a par- 
ticular configuration of impurities. A typical Feymann 
diagram contributing to the above expression, is a "bub- 
ble" diagramm with an electron line going from r to r' 
over a given set of impurities, and then back from r' to 
r oven a (generally) different set of impurities; this con- 
tribution to Eq. (|^) contains a phase factor 



exp{ifc(z) [L(r ^ r') + L(r' ^ r)]} , 



(2) 



where k{z) is the (complex) wavevector corresponding to 
the complex energy z, and where L{r — > r') (L(r' — r)) 
is the length of the path from r to r' (r' to r) over the 
corresponding set of impurities. 

For points r and r' not very close to each other, when 
summing over all diagrams (i.e., over all paths over im- 
purities), the above factor oscillates rapidly, leading to a 
strong cancella|tion. Thus, as pointed out by Bulaevskii 
and Panyukov,0 the only significant terms are due to the 
paths going through impurities lying on a straight line be- 
tween r and r' (the impurities being passed sequentially, 
without back tracking); for all such paths, the phase fac- 
tor is 



exp[2ifc(z)|r- r'l 



(3) 



We wish to stress that the present situation is com- 
pletely different from the one encountered when comput- 
ing the two-point conductivity cr(r, r'), which involves the 



product G(r, r'; £f + iO+)G(r', r; ep - iO+). This yields 
an oscillatory factor 



exp {ikp [L{r ^ r') - L{r' r)]} . 



(4) 



Because of the minus sign in the above expression, there 
is a significant contribution from all diagrams such that 
the path r' ^ r is the reverse of the path r — > r', so 
that L{r — > r') = L{r' — > r); these are the ladder dia- 
grms characteristic of a diffusive process, which give the 
leading contribution to cr(r,r'). j-, 

Following Bulaevskii and PanyukovJ3 we obtain the ex- 
pression of the coupling: 



J(r,r') 



cos[2A:j;-|r — r' 



(r,r')] 



r — r 
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the phase shift due to impurity scattering is 

-2 /■I'-'"'! 
(/)(r,r') = - — / dsU{r + ns), 
hvp Jo 



(5) 



(6) 



where n is a unit vector of the (r, r') axis and U{r) is the 
impurity scattering potential. Equations (^) and (||) are 
valid if the perturbation potential U{r) is small compared 
to the Fermi energy. As there is a specific distribution 
of impurities between a pair of ions there is a definite 
phase; the coupling is phase shifted but it is not damped 
by impurity scattering (its decay law is as |r — r'|~'^ like 
in the pure system). 

The distribution of J(r,r') is thus determined by the 
distribution of the phases (j>{r,r'). To ascertain the dis- 
tribution of phase angles we fill the space between two 
spins with cubes of length a, and rewrite the phase inte- 
gral Eq. (^ as a sum over the N cells {N = \r — r'\/a) 
crossed by the trajectory 



N 

1=1 



(7) 



'2Ua/hvF (if impurity in cell i) 

(if no impurity in cell i) 



If an impurity is in a cell it yields a contribution k = 
—2Ua/hvp — —{Sny/'^U/Ep to the phase; otherwise it 
gives zero. For a given concentration of impurities c there 
is a probability c at each site of their being an impurity 
and therefore of picking up a phase of k, and a probability 
of 1 — c of picking up zero. Therefore the phase in Eq. (|^) 
is just a binomial distribution of N events. 

^From the distribution of phases, one can compute the 
various moments of the exchange interaction distribution, 



J"(i?)^(J"(r,r')), 



(8) 



where the angular brackets indicate that we average over 
all possible configurations of impurities. The averag- 
ing restores the translational and rotational invariance, 
therefore the moments J" depend only on the distance 
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R = \r — r'l betweeii|-the two spins. As was emphasized 
by various authors,Bu the first moment J{R) bears Uttle 
physical significance in the case of magnetic ions embed- 
ded in a disordered nonmagnetic host; for instance, the 
transition temperature of spin glasses is determined by 
the second moment J^(i?), not by the first moment. To 
illustrate this point, we compute now the first two mo- 
ments, J{R) and J'^{R). 

The configuration averaged exchange interaction (first 
moment) is given by 



Re 



JiR) 



^2ikpR ^gi,/,(r,r')^j 



(9) 



The characteristic function the phase,0 for the bino- 
mial distribution is given aa 



[(1 



I R/a 



(10) 



where we have used R — Na. 

By taking the logarithm of the characteristic function 
for the binomial distribution we find in the limit of low 
impurity concentrations 



,i0(r,r') 



exp [i 



c e 



1) R/a] 



c(s'mK,)R/a —c{l~cos K)R/a 



(11) 



Thus, we obtain an exponential decay of the exchange 
interaction with a decay length A given by 



A = - (1 — cos k) « ^ 
a 



V3 



2a \eF 
and a shift in the Fermi wavevector 



Sk ^ 



c sm K 
2a 



-c (37r^)'^' 
2a 



U 

SF 



(12) 



(13) 



The wavevector shift Sk is simply due to the shift in the 
average value of the potential; in the following it will be 
incorporated into a redefinition of the Fermi wavevector 
kp- 

In this way we find that the average RKKY coupling 

is 



JiR) 



cos(2fci?i?) e 

i?3 



(14) 



This result was first obtained by de Gennes.0 

The second moment. ]s. easily calculated directly from 
Eq. (^ and we obtairo" 



JHR) 



1 1 

2 R^' 



Thus it appears clearly that 

j^{r)-j\r) 
1\r) 



(15) 



(16) 



i.e., that the exchange interaction between a pair of mag- 
netic ions is not a seitaveraging quantity (in the sense of 
Kohn and Luttingert2l) . 

In the case of a system with a non-spherical Fermi 
surface, the configuration averaged exchange interaction 
takes a form similar to Eq. (|lj), but the wavevector of os- 
cillations and the decay length both depend on the direc- 
tion n. We stress that the decay length of the exchange 
interaction corresponding to a particular direction gener- 
ally has no simple relation with the transport mean free 
path, because the later results from averaging over all 
directions. 



III. TWO SHEETS OF SPINS 
A. Perfectly flat layers 

We now consider the interlayer exchange coupling be- 
tween two ferromagnetics layers, -Fi and F2. These are 
modelled by taking two sheets of magnetic ions of nor- 
mal coordinates rj_i and rj_2, respectively. Within a 
given sheet, we assume that all the magnetic moments 
are maintained parallel to each other by some intralayer 
exchange coupling (which we do not describe explicitely 
here); thus, the only variable is the angle between the 
magnetizations of the two sheets. 

Following Yafet|lj we express the coupling between 
Fi and ^2, as the sum over the pairs of magnetics ions 
(ri, r2) (divided by the total area 5), with ri = (ryi, r_Li) 
belonging to Fx and V2 = (r||2,r_L2) belonging to F^. 
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d^r||i / d^rp J(ri,r2). 



(17) 



To compute this, we first sum over all pairs (ri, r2) with 
ri — r2 parallel to a given direction, i.e., we rewrite the 
above equation as 



Iir±i,r±2) = Jd^P\\ K{r±i,r±2;P\\), 



(18) 



with 



K{rj_i,r±2\P\\) ^ ^ j d?Y\\ J (^(r,, , r^i), (r,| p,, , r^a^ 

(19) 

It is easy to see that, when summing over r||, all con- 
figurations of impurities between ri = (r||,7-j^i) and 
^2 = (r|| -l- P||,r_L2) are encountered, thus this is equiva- 
lent to performing a configuration average, i.e.. 



K{r^i,r±_2;p^) = J{Ri,2), 



where 



i?i,2 = JpI + Z?2 



(20) 



(21) 
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is the distance between the spins for the pairs considered, 
and D = |rj_i — r±^2\ is the distance between the two 
sheets of spins. 

This impHes that, in contrast to the exchange interac- 
tion between two magnetic ions, the exchange coupling 
between two sheets of spins is self-averaging in the sense 
of Kohn and Luttinger.t^l Thus, it is only a function of 
the distance D between the two sheets, and not of rj^i 
and r_L2 separately. 

The remaining integration over p|| is then easily calcu- 
lated by using Yafet's method^ and we find 



I{D) - 27rRe 



7 



with 7 = 2fci? + iA ^; finally, wc obtain 

sin (2fcf £) + ip) rj/^ 



with 



Lp = arctan 



/ -1 



-1 



(22) 



(23) 



(24) 



Thus, the presence of impurities in the non-magnetic 
spacer layer leads to an exponential decay of the inter- 
layer exchange coupling with the distance between mag- 
netic layers; this result is in contrast to the one obtained 
for the exchange interaction between magnetic ions in the 
previous section. It is traced back to the self-averaging 
character of the interaction between planes. Another ef- 
fect of impurity scattering is the phase shift p [Eq. (p^j . 

The results of the present subsection are in full agree- 
ment with the ones obtained previously by Bruno et al. 
by using first-principles calculart|ktfis together with the 
"vertex cancellation theorem. "tSEj Moreover, this pre- 
vious study allows us to generalize the above result to 
the case of a spacer material with non-spherical Fermi 
surface; in this case, one obtains 



(25) 



In the above equation, q±a and are the real and 
imaginary parts of sta|tipnary spanning vectors of the 
complex Fermi surfac£3 of the alloy spacer material 
(since there may be several such vectors, they are la- 
beled by the index a); la and 0^ are the corresponding 
amplitude and phase. 

The model calculation presented here provides a sim- 
ple physical explaination for the "vertex cancellation the- 
orem." As explained at the beginning of Sec.||, only the 
paths going in straight line between two ions contribute 
significantly to the exchange interaction between them; 
this, together with the self-averaging property for the 
coupling between layers, forms the physical basis of the 
"vertex cancellation theorem." 



B. Effect of interface roughness 

Next we discuss the effect of interface roughness on the 
interlayer exchange coupling. To be specific, we consider 
the case where the normal coordinates r±i and r±2 char- 
acterizing Fi and F2 vary with the in-plane coordinates 
rp and r||2. The roughness is characterized by at least 
two parameters: the average amplitude of the fluctua- 
tions of r±i and r±2, and the lateral correlation length. 

The simplest approach to the effect of roughness con- 
sists in calculating the effective interlayer exchax|gc_cou- 
pling by averaging over thickness fluctuations .□'EZfll3 In 
order to be allowed to do so however, some conditions 
must be satisfied. The first condition is that the lateral 
correlation length of the roughness, ^, should be large 
enough for the interlayer exchange coupling to be kt 
cally well defined; typically this requires that ^ > Z?.E£l 
On the other hand, we wish to consider that the sheets 
of spin are uniformly magnetized; but local fiuctuations 
of spacer layer thickness induce local interlayer coupling 
fluctuations, that tend to produce local fluctuations of 
the magnetization direction in the magnetic layers. Thus, 
in order to keep the magnetization direction constant 
in the magnetic layers, the intralayer exchange coupling 
must be large enough, and the correlation length ^ small 
enough; in practice, this condition ^ is not very restric- 
tive, and we shall not consider it further. 

Within the above conditions, the effective interlayer 
exchange coupling / is given by 



dD P{D) I{D), 



(26) 



where P[Fi)-. is the distribution function of spacer 
thicknesses. Thus, we have 



/ - 27rRe 



I dD P{D) 



eirD 



(27) 



if the width of the distribution of thicknesses is small 
compared to the average thickness D, then the above 
equation becomes 



27rRe 



-i e'''^ 
- ^(7) — 
7 D 



where 



A(7) = [dD P{D) e'''(^-^) 



(28) 



(29) 



is the form factor for the roughness. For a Gaussian 
distribution of width ct, i.e.. 



P{D) = 



1 



,-{D-dY/2c 



2TTa 



(30) 



one has 
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-2kf 



Thus, the effective couphng becomes 
T ^ „-2fcFV^ sin (2^^!? + V3 + V) 



/A 



where 



A 



(31) 



(32) 



(33) 



The effect of interface roughness is essentially to attenu- 
ate the oscillatory coupling by the factor exp (— 2A:|,(t^). 
In contrast to the attenuation due to impurities in the 
spacer layer, this effect is independent of the thickness 
of the spacer layer, D. The phase shift ip in Eq. (^) 
is a combined effect of the roughness and of the impu- 
rities. An alternative approach for treating the effect of 
roughness is also presented in the Appendix. 

In the more general case where the interlayer exchange 
coupling comprises several oscillatory components, the 
effect of the roughness is to strongly attenuate those os- 
cillatory components which have a period of the |Qp4(?f, of 
(or smaller than) the amplitude of roughness, cr.l 



IV. CONCLUSION 

By using a semi-classical approach, we have shown how 
the various kinds of defects, i.e., impurity scattering in 
the bulk of the layers and interface roughness, modify the 
exchange coupling between a pair of magnetic ions, or be- 
tween two magnetic sheets. Our approach emphasizes the 
intrinsic self-averaging character of the coupling between 
layers, by contrast to the interaction between ions. 

We stress that although the same type of defects (im- 
purities, roughness) play an essential role both in the 
interlayer exchange coupling and in the transport prop- 
erties of magnetic multilayers, one cannot draw any sim- 
ple relationship between their influence in the two cases. 
For example, in the absence of imBurity scattering, it has 
been shown by Zhang and LevyB that the effect of in- 
terface roughness on the perpendicular transport can be 
described in terms of an effective mean free path; but, as 
the above discussion has shown, the interlayer exchange 
coupling is not exponentially damped by the roughness 
alone, and it would be wrong to believe that this effec- 
tive mean free path is of any significance for the interlayer 
exchange coupling. 

Here, we wish to comment on the limitations and ex- 
tensions of the results obtained in this paper. For the 
sake of clarity and simplicity, we have restricted our- 
selves here to a very simple model (magnetic layers of 
infinitesimal thickness), and to a perturbative approach 
(expressing the exchange interaction via a susceptibil- 
ity). In order to treat more realistic systems, with mag- 
netic layers of finite (or infinite) thickness, a more sophis-, 
ticated approach, such the one developed by Bruno,ta 



should be used; however, we expect that the conclusions 
obtained here would still hold. The effect of impurity 
scattering is described in tewns-|Of complex wavevectors 
and complex Fermi surfaceE3'cB Impurity scattering in 
the spacer layer gives rise to an exponential damping of 
the coupling, but not impurity scattering in the magnetic 
layers (in this case, only the amplitudes and phases are 
affected) , which agrees with what one would expect intu- 
itively. Interface roughness modifies not only the spacer 
layer thickness, but also the thickness of the magnetic 
layers. It is known that the interlayer exchange coupling 
varies not only with the spacer thickness, but also with 
the thickness of the magnetic layer;E£l however, the latter 
is a secondary effect and can be neglected here. Thus, we 
expect that the effect of roughness would be essentially 
the same as described within the simple approach of the 
present paper. 

Finally, we have considered only the limit case of "ge- 
ometrical" roughness. An opposite limit case is the one 
of an interdiffusion of the magnetic material and spacer 
material near the interface. In such a situation, it is com- 
pletey inappropriate to discuss the effect of roughness in 
terms of fluctuations of the-apacer thickness. Rather, as 
done by Kudrnovsky et aL,cJ one can consider that there 
is a thin layer of disordered magnetic alloy in the inter- 
face region, whose magnetization remains parallel to the 
magnetization of the magnetic layer nearby. In such a 
case, as one would expect intuitively, one finds that the 
interdiffusion modifies the amplitude(s) and phase(s) of 
oscillatory coupling, but does not alter the period(s) nor 
the decay. 
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APPENDIX: 



As we have shown in Sec. IIIB interface roughness 
makes it necessary for us to average the interlayer cou- 
pling over a distribution of spacer thicknesses; see Eq. |2^. 
Here, we present an alternate way of arriving at the at- 
tenuation factor due to roughness (thickness fluctuations) 
by using the canonical transformation introduced by 
Tesanovic, Jaric and Maekawa to transform a film with a 
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rough boundary to one with smooth ones,iil and extended 
to films witla-|two rough boundaries by Meyerovich and 
Stepaniants.ca In doing this the transformation induces 
a perturbation (scattering potential) into an otherwise 
impurity free layer. It follows that we can replace the 
average over spacer layer thicknesses by one over a fluc- 
tuating phase shift that is induced by the canonical trans- 
formation that replaces the spacer with rough boundaries 
(interface roughness) by one with smooth boundaries. 

By following Refs. ^,^or Trivedi and AshcroftcJ the 
perturbation due to the surface roughness is 



Kiurfaco(r) = ^?/(r||) {[f_LP_L ^0 



c.c.} 



^2ry(r||)i7o 



where 



2m 2m 



D '■ 



(Al) 

(A2) 
(A3) 



and Vo{r±) is the confining potential, SD{r\\) is the vari- 
ation in thickness over the surface, and D is the average 
thickness of the layer. 

As our treatment up till now has neglected the con- 
fining potential, or equivalently in the limit of large D, 
the distribution of eigenvalues is quasi-continuous. The 
matrix element of the perturbation, Eq. (Al), between 
states of Hq at the Fermi level is 



(k|| , n I KiurfaccI k|| + q[| , n) = 2?7(q|| )eF 



(A4) 



where n labels the states referring to energy levels in the 
r± direction, and f)(q||) is the two-dimensional Fourier 
transform of ?y(r||). 

By placing this scattering potential in the semiclassi- 
cal expression for the Green's functional, we find the ex- 
change interaction between two spins located on Fi and 
F2 is given again by Eq. (^) , but with a fiuctuating phase 
shift, Eq. (i), of 



(ri, rz) = 2kFRi,2 '7(r!|i - r||2)- 



(A5) 



By proceeding as in Sees. and Q we obtain an expres- 
sion of the interlayer coupling averaged over thickness 
fluctuations: 



/ - 27r Re 



_• i2kFD 

A'{2kp) 
2kF D 



(A6) 



those obtained in Sec. [I], when no impurities are present, 
i.e., = 0. 
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A'{2kF) 



^ /i2kFSD{v^^) 



(AT) 



As in Eq. ( [l3| ) we include a shift Sk by redeflning fc^. 
Comparison with Eq. ( |29| ) shows that A'{2kF) is equal to 
yl(2fci?). Thus, one obtains results that are equivalent to 
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